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■^ , Abstract: We consider the zero dissipation limit of the full compressible Navier-Stokes 

,^ ' equations with Riemann initial data in the case of superposition of two rarefaction waves and a 

j^ ■ contact discontinuity. It is proved that for any suitably small viscosity e and heat conductivity 

K satisfying the relation ()1.3p . there exists a unique global piecewise smooth solution to the 
compressible Navier-Stokes equations. Moreover, as the viscosity £ tends to zero, the Navier- 
Stokes solution converges uniformly to the Riemann solution of superposition of two rarefaction 
^ . waves and a contact discontinuity to the corresponding Euler equations with the same Riemann 

^^ I initial data away from the initial line i = and the contact discontinuity located at x = 0. 

(N 

^ 1 Introduction 

o 

CN I We study the zero dissipation limit of the solution to the Navier-Stokes equations of a 

compressible heat-conducting gas in Lagrangian coordinate: 



Vt-Ux = 0, 

Uxs 

~ (1.1) 



Ut+Px = e( — )x, 
2 ^ 

e + Tr)i + (P'^h = k{—)x + e( )x 

2 V V 



with Riemann initial data 
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where the functions v{x, t) > 0, u{x, t), 9{x, t) > represent the specific vohime, velocity and the 
absolute temperature of the gas, respectively. And p = p{v, 0) is the pressure, e = e{v, 6) is the 
internal energy, e > is the viscosity constant and k > is the coefficient of heat conduction. 
Here we consider an ideal and polytropic gas, that is 

RO Re 

P 



V 7 — 1 

with 7 > 1, i? > being gas constants. 

The study of the asymptotic behavior of viscous flows, as the viscosity tends to zero, is one 
of the important problems in the theory of compressible fluid flows. When the solution of the 
inviscid flow is smooth, the zero dissipation limit problem can be solved by classical scaling 
method. However, the inviscid compressible flow contains discontinuities, such as shock waves, 
in general. In this case, it is also conjectured that a general weak entropy solution to the inviscid 
flow should be the strong limit of the solution to the corresponding viscous flows with the same 
initial data as the viscosity vanishes. 

It is well known that the solution to the Riemann problem for the Euler equations consists of 
three basic wave patterns, that is, shock, rarefaction wave and contact discontinuity. Moreover, 
the Riemann solution is essential in the theory for the Euler equations as it captures both local 
and global behavior of general solutions. 

For hyperbolic conservation laws with the uniform viscosity 

ut + f{u)x = euxx, 

where f{u) satisfies some assumptions to ensure the hyperbolic nature of the corresponding 
inviscid system, Goodman-Xin [4] verified the limit for piecewise smooth solutions separated 
by non-interacting shock waves using a matched asymptotic expansion method. Later, Yu |33j 
proved it for hyperbolic conservation laws with both shock and initial layers. In 2005, important 
progress made by Bianchini-Bressan[T] justifies the vanishing viscosity limit in BV-space even 
though the problem is still unsolved for the physical system such as the compressible Navier- 
Stokes equations. 

For the compressible isentropic Navier-Stokes equations where the conservation of energy in 
p.ip is neglected in the isentropic regime, Hoff-Liu [TT| firstly proved the vanishing viscosity 
limit for a piecewise constant shock with initial layer. Later, Xin [31] justified the limit for 
rarefaction waves. Then, Wang [29] generalized the result of Goodman-Xin [3] to the isentropic 
Navier-Stokes equations. 

Recently, Chen-Perepelitsa [2] proved the convergence of the isentropic compressible Navier- 
Stokes equations to the compressible Euler equations as the viscosity vanishes in Eulerian coor- 
dinates for general initial data by using compensated compactness method if the far field does 
not contain vacuum. Note that this result allows the initial data containing vacuum in the inte- 
rior domain. However, the framework of compensated compactness is basically limited to 2 x 2 
systems so far, so that this result could not be applied to the full compressible Navier-Stokes 
equations (jl.ip . 

For the full compressible Navier-Stokes equations, there are investigations on the limits to 
the Euler system for the basic wave patterns in the literature. We refer to Jiang-Ni-Sun [17] and 



Xin-Zeng [32] for the rarefaction wave, Wang [30] for the shock wave, Ma [21] for the contact 
discontinuity and Huang- Wang- Yang [141 [TS] for the superposition of two rarefaction waves and 
a contact discontinuity and the superposition of rarefaction and shock waves, respectively. We 
should point out that the limit shown in [17] was for the discontinuous initial data while the 
other results mentioned were for (well-prepared) smooth data. 

In this paper, we shall investigate the zero dissipation limit of the full Navier-Stokes equations 
(jl.ip with Riemann initial data ()1.2p in the case of the superposition of two rarefaction waves 
and a contact discontinuity. The local and global well-posedness of the full system (jl.ip or the 
corresponding isentropic system with discontinuous initial data is systematically studied by Hoff, 
etc. , see [3 [6l [3 [HI [H [lOl [3] . In order to get the zero dissipation limit to the Riemann solution 
of the Euler system, we shall combine the local existence of solutions with discontinuous data 
from [7| and the time-asymptotic stability analysis to the compressible Navier-Stokes equations 
()2.2p . Compared with the previous result [H] where the same limit process is studied for (well- 
prepared) smooth initial data, the main difficulty in the proof here lies in the discontinuity of 
the initial data. The discontinuity of the initial data for the volume v(t, x) will propagate for all 
the time along the particle path due to the hyperbolic regime while the smoothing effects will 
also be performed on the velocity u{t, x) and the temperature 6{t, x) by the parabolic structure, 
and this interaction of the discontinuity and smoothing effects brings technical difficulties. To 
circumvent such difficulties, we shall choose suitable weight functions to carry out the weighted 
energy estimates in terms of the superposition wave structure (see Remark 3.7), and use the 
energy method of Huang-Li-Matsumura [12] for the stability of two rarefaction waves with a 
contact discontinuity in the middle, where the authors obtained a new estimate on the heat 
kernel which can be applied to the study of the stability of the viscous contact wave in the 
framework of the rarefaction wave (see Lemma 13. 6p . Namely, the anti-derivative variable of the 
perturbation is not necessary and the estimates to the perturbation itself are also available to 
get the stability of the viscous contact wave. 

Without loss of generality, we assume the following relation between the viscosity constant 
e and the heat-conducing coefficient k, of system (jl.ip as in [T7| : 

K = 0{e), as e — > 0; 

u = > c > for some positive constant c, as e — ?• 0. 

e 

If K = e = in (|l.ip . then the corresponding Euler system reads as 

'' Vt-Ux=0, 

u \ 

It can be easily computed that the eigenvalues of the Jacobi matrix of the flux function to (jl.4p 
are 

Ai = -W — , A2 = 0, A3 = W — . (1.5) 

\ V \ V 

It is well known that the flrst and third characteristic fields of ()1.4p are genuinely nonlinear and 
the second one is linearly degenerate (see [28]). 



For the Euler equations, we know that there are three basic wave patterns, shock, rarefaction 
wave and contact discontinuity. And the Riemann solution to the Euler equations has a basic 
wave pattern consisting the superposition of these three waves with the contact discontinuity in 
the middle. For later use, let us firstly recall the wave curves for the two types of basic waves 
studied in this paper. 

Given the right end state (v+,n+,0-)_) with f+,^+ > 0, the following wave curves in the 
phase space {{v,u,6)\v > 0,9 > 0} are defined for the Euler equations. 

• Contact discontinuity curve: 

CD{v^, u+,^+) = {{v,u,6)\u = Uj^,p = p+,v ^ «+}. (1-6) 

• i-Rarefaction wave curve {i = 1,3): 



Ri{v+,u+,e+) := { {v,u,e) 



u<u+,u = u+- Xi{r],s+)dr], s{v,6) = s+>, (1.7) 



V ^ 




/ w^i + v"'^ + v''^ 


u 


{t,x) = 


^n _^ ^cd _^ ^ra 


e > 




V 0^1 + 9'"^ + 0^3 



where s+ = s(u+,6'+) and A^ = Xi{v,s) defined in (jl.Sp is the i-th characteristic speed of the 
Euler system (II. 4p . 

Now, we define the solution profile that consists of the superposition of two rarefaction 
waves and a contact discontinuity. Let {v-,U-,6-) £ Ri-CD-R3{v+,u^,9^). Then, there exist 
uniquely two intermediate states ( v* , u^, , 6*=^ ) and {v*,u*,9*), such that {v^:,u^,6^) £ Ri{v-,u-,9-) 
{v^,u^,9^) £ CD{v*,u*,9*) and_{v*,u*,9*) £ R'i{v+,u+,9+). 

Thus, the wave pattern {V , U, Q){t, x) consisting of 1-rarefaction wave, 2-contact discontinu- 
ity and 3-rarefaction wave that solves the corresponding Riemann problem of the Euler system 
(jl.4p can be defined by 

'"■--''"- ^ / v, + v* \ 

{t,x)- u^+u* , (1.8) 

V ^* + ^* / 

where {v^^ , u'^^ , 9^^){t, x) is the 1-rarefaction wave defined in (|1.7p with the right state (f+, u_|_, ^_|_) 
replaced by (f =„ , n^, , 0* ) , {v'^'^,u'^'^,9'^'^){t,x) is the contact discontinuity defined in ()1.6I) with the 
states (f_,n„,0_) and (^+,^+,0+) replaced by (z;* , n^, , 6** ) and {v*,u*,9*) respectively, and 
{v^^,u'^^,9^''^){t,x) is the 3-rarefaction wave defined in (II. 7p with the left state {v-,U-,9-) re- 
placed by {v*,u*,9*). 

Now we state the main result as follows. 

Theorem 1.1. Given a Riemann solution {V, IJ, ©)(t, x) defined in (jl.Sp . which is superposition 
of two rarefaction waves and a contact discontinuity for the Euler system ()1.4p . there exist small 
positive constants 5o and Eq, such that if e < Sq and the wave strength 5 = \{v+—V-, u+—u-,9+ — 
d~)\ ^ ^0, then the compressible Navier-Stokes equations (II. ip with (II. 2p and (II. 3p admits a 
unique global piece-wise smooth solution {v^,u^,9'^){t,x) satisfying that 

• The quantities u^,9^, p{v^,9^) — e^ and :^ are continuous for t > 0, and the jumps in 
v^ ,u%,9% at X = satisfies 

|(K(t,O)],[<.(t,O)],[0^.(t,O)])|<Ce-f, 
where the constants C and c are independent of t and e. 



• Moreover, under the condition (jl.3p . it holds that 

lim sup \{v^,u^,e^){t,x)-{V,U,e){t,x)\=0, V/i > 0, (1.9) 

where S/, = {{t,x)\t > /i, -^ > /ie°,0 <a< i}. 

Remark 1.2. Theorem M . 1\ shows that, away from the initial time t = and the contact discon- 
tinuity located at x = 0, there exists a unique global solution (w",ii^,^^)(t, x) of the compressible 
Navier-Stokes equations (jl.ip which converges to the Riemann solution {V , U, Q){t,x) consisting 
of two rarefaction waves and a contact discontinuity when e and k satisfy the relation (|1.3p and 
£ tends to zero. Moreover, the convergence is uniform on the set S/j for any h > 0. 

Notations. In the paper, we always use the notation t = / + / i II ' II to denote the 

... .,K, „o..n, ,M .„ de„o.e .e .ecewise .^ J: .^ Hrfl I ;.,, , , .. 

Jn 

II • jji represent the H^(R) norm and piece- wise i?^(R ) norm, respectively. And the notation 

[•] represents the jump of the function • at x = or y = if without confusion. 

2 Approximate profiles 

Introduce the following scaled variables 

y = ^, r = ^, (2.1) 

and set 

(^;^n^r)(t,x) = (^;,n,e)(T,y). 

Then the new unknown functions {v,u,6){T,y) satisfies the system 

Vr — Uy = 0, 

Ur+Py-[ — )v, ^2.2) 

-Or + pUy = l^{ — )y + — , 
1 V V 
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with the scaled heat conductivity i^ = - in (|1.3p satisfying 

1^0 ^ ^ ^ ^1) uniformly in e as e — )■ 0+, for some positive constants z/q and vi. 

Note that the Riemann solution {V, U, @){t, x) in (jl.Sp is invariant under the scaling trans- 
formation ()2.ip . thus to prove the limit ()1.9p in Theorem ll.il it is sufficient to show the following 
limit 

hm sup \{v,u,9){T,y)-{V,U,e){T,y)\=0, V/i > 0, (2.3) 



where Sj^ is the corresponding region of S/^ in the new coordinates (r, y) defined by 



4 = {{T,y)\T > - -y^ > ^,0 < a < 1). (2.4) 



Now we study the Navier-Stokes equations (I2.2p . The corresponding wave profiles to (I1.6P 
and (jl.7p can be defined approximately as follows. We start from the viscous contact wave to 
([LSD. 

2.1 Viscous contact wave 

If (f„,u_,0_) € CL'('y+, n+, 0+), i.e., 

u- = u+, p- = p+, V- ^ v+, 
then the Riemann problem, that is, the Euler system (II. 4p with Riemann initial data 

(t;_,n_,6'_), y < 0, 



{v,u,e){T = ^,y) - 
admits a single contact discontinuity solution 

{v^'^X',e^'){r,y) = I 5^-'^+'^-)' 2/ < n+r, r > 0, 

^ ^^ ^ \ (f+,n+,0+), y>u+T,T>0. 

As in [13], the viscous version of the above contact discontinuity, called viscous contact wave 
{V^^ , U^^ , 0'^^)(r, y), can be defined as follows. Since it is expected that 

P'^^ ^p+=p^, and |[/'^^-n+| < 1, 

the leading order of the energy equation (]2.2P 3 is 

7? f) 

-jer+P+Uy = lj{-^)y. 



7 
Then, similar to |12j or [T3], one can get the following nonlinear diffusion equation 



e. = a -^ , e(r,±) = 0±, a 



The above diffusion equation has a unique self-similar solution G(r, y) = 0( ,y_ ). 
Thus, the viscous contact wave {V'^'^, U^^, @^^){T,y) can be defined by 

p+ 

f/<^-(r,!,)=., + ;%^|^, (2.5) 

e°°(r.rf = e(r.,) + ^^-'''-'-'> e.. 

7P+ 



Here, it is straightforward to check that the viscous contact wave defined in (j2.5p satisfies 



e - 0±\ + {1 + T)2\Qy\ + {1 + T)\eyy\ =0(1)5^^6"^+^, as \y\^+oo, 



(2.6) 



where 6 = |0+ — 0_| represents the strength of the viscous contact wave and cq is a positive 
constant. Note that in (|2.5p . the higher order term "')'~'^V7- ) q^ jg introduced in G'^^(r, y) to 
make the viscous contact wave {V'-'^ , U^^ , 0'^^)(t, y) satisfy the momentum equation exactly. 
Correspondingly, {V^^ ,U ,Q'~^^){T,y) satisfies the system 



where P'^^ 





ttCD , uCD f-^y \ 


« 


W 1 ^y yyCDJy^ 

p ,ptCD^ (TtCD\2 
^ (P)CD 1 pCDrrCD _ (^y \ . ^'^y > 


RQ^ 


I 7-1^- ' ^ '^y 'WcDjy 1 VCD 
— and thfi prrnr tfirm QCD Rafisfifis 



+ Q 



CD 



yCD 



(2.7) 



Q 



CD 



0(1)5^^(1 + r)-^e-^+^ 



as jyj 



+00, 



for some positive constant cq. 



(2.8) 



2.2 Approximate rarefaction waves 

We now turn to the approximate rarefaction waves to ()1.7p . Since there is no exact rarefaction 
wave profile for the Navier-Stokes equations, the following approximate rarefaction wave profile, 
which satisfies the Euler equations, is motivated by [31]. For the completeness of presentation, 
we include its definition and the properties in this subsection. 

If (?;_ , u_ , 0_ ) G i?j(t'_|_, n+, 0+), (i = 1, 3), then there exists an i-rarefaction wave (f^%«^S 0^*)(y/T) 
which is a global solution of the following Riemann problem: 



u,, 



0, 



Ur +Py{v,9) = 0, 
n 

-Or +p{v,6)uy = 0, 

(v_,u_,6'_), y<0, 
{v+,u+,9+), y > 0. 
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{v,u,0){O,y) 



Consider the following inviscid Burgers equation with Riemann data: 



(2.9) 



Wr + WW^ 



y 



w{t = 0,y) 



W-, y<0, 
w+, y>0. 



(2.10) 



If W- < W-^- , then the Riemann problem (j2.10p admits a rarefaction wave solution 



T 






(2.11) 



Thus, the Riemann solution in (|2.9p can be expressed explicitly through the above rarefaction 
wave (|2.1ip to the Burgers equation, that is. 



' /■^(T,y) = s{v''^{T,y),9''^{T,y)) = s+, 

w± = Xi± := Xiiv±,e±), 

«;'-(^) = A,K'(r,y),s+), 
r 

u'''{T,y) = u+- Xi{v,s+)dv. 



(2.12) 



In order to construct the approximate rarefaction wave (y^% C/^% 0^')(r, y) corresponding 
to (|1.7p . we first consider the following approximate rarefaction wave to the Burgers equation: 



Wr + WWy = 



y — "' 



, W+ + W- Wj^ — W- 

[ w{0,y) =wo{y) = \ tanhy. 

Note that the solution w {T,y) of the problem (12.13P is given by 

w^{t, y) = wq{xq{t, y)), x = xo(t, y) + wo{xo{t, y))T. 
And w^{t, y) has the following properties, the proof of which can be found in \2,2\ [31]: 
Lemma 2.1. Let w^ < W-^, then ()2.13p has a unique smooth solution w {T,y) satisfying 

(1) w. < w^{T,y) < w+, {w^)y{T,y) > 0; 

(2) For any 1 < p < +00, there exists a constant C such that 

\\-l^wR{T,-)\\Lnn)<Cm\n{{w+-w.), {w+ - w.YIpt~^+^Ip], 



(2.13) 



4.«^^(T,-)llLnR)<Cmin{( 



w+ — w. 



0, r-^}; 



(3) Ify-W-T<0, then 



\w^{T,y)-W-\ < {w+-W-)e-'^\y-'"-''\, 
\^w''{T,y)\<2iw+-w.)e-^\y-^--\- 



If y — w+T > 0, then 



-2\y-w+T\ 



\w^{T,y) — w^\ < (tt;+ — W-)e 
li^^(r,2/)| < 2{w+ - w-)e-^\y-^+-\- 



(4) sup \w^{T,y) — w^{^)\ < mill {tt;+ — W-, ^ ln(l + r)|. 
yen 

Then, corresponding to ()2.12p , the approximate rarefaction wave profile denoted by (1/ ' , [/ ' , * ) (r, y) {i 
1,3) to (|1.7|) can be defined by 

w± = Ai± := h{v±,0±), 

«;^(l + r,y) = \(y^»(r,y),SH.), (2.14) 

U^'{T,y) = u+- Xi{v,s+)dv. 



Note that (F^S C/^S 9-^')(r, y) defined above satisfies 



^ y/^' - C//^» = 0, 



U^' + P^' = 0, 



R 



(2.15) 



-ef' + p^'C//^' = 0, 



I 7-1 
where P^' = p{V^%Q^'). 

By virtue of Lemmas 12 .11 the properties on the approximate rarefaction waves iV "^ , C/ ' , * ) (r, y) 
can be summarized as follows. 

Lemma 2.2. T/ie approximate rarefaction waves (y^\ C/^*, 0^*)(t, y) (z = 1,3) constructed in 
()2.14p /lave the following properties: 

(1) U^^{T,y)>^fory^n,T>Q; 

(2) For any 1 <p < +oo, the following estimates holds, 

||(y^',^^',G«Oj/llLp(d,) <Cmin{<5«% (^^^'/^l + ^)~'+'/^}, 
||(y«',t/«',e«Oj/yllLP(d,) <Cmin{<5«% (l + r)"!}, 

where 5 * = |(v+, v_,u+,u_,0+, 0_)| zs i/ie i-rarefaction wave strength and the positive 
constant C is independent of r, but may only depend on p and the wave strength; 

(3) Ify>\i+{l + T), then 

I (y^i , f/^i , e^i ) (r, y) - (t;„ , n_ , 0_ ) I < C5^i e-2|s'-^i+ (1+^)1 , 
\{V^\U^\Q^'%{T,y)\ < c;5fiie-%-^i+(i+^)l; 
//y< A3-(l + r), then 

\{V^^ ,U^^ ,Q^-^){T,y) - {v+,u+,e+)\ <c;5R3e-2|j/-A3-(i+r)|^ 

1(1/^3 ^f/fi3^ 9^3 )^(^^y) I < C'5i?3e-2|y-A3-{l+r)|. 

(^^j There exists a positive constant C , such that for all r > 0, 



>y^ 



C 



sup|(y^%c/^%e^O(^,y)-(^''%^''%^''0(-)l < in(i + r). 

2/eR T 1 + T 



2.3 Superposition of rarefaction waves and contact discontinuity 

Corresponding to (jl.Sp . the approximate wave pattern (V,U,Q){T,y) of the compressible 
Navier-Stokes equations (12. 2p can be defined by 




yRi _^ yCD _^ yRs 

(r, y) = ( [7^1+ U^^ + C/^3 

qRi + eCD ^ 0R3 



(^>y) 




(2.16) 



where {V^^ ,U^^ ,Q^^){T,y) is the approximate 1-rarefaction wave defined in (j2.14p with the 
right state (v+,u+,0+) replaced by (i;* , tt* , 0* ) , {V ,U ,Q ){T,y) is the viscous contact 
wave defined in (j2.5p with the states (t;_,u_,0_) and {v+,u+,6+) replaced by {v^,u*,6^) and 
{v* ,u* , 0*) respectively, and iV -^ , U -^ , & '^){t, y) is the approximate 3-rarefaction wave defined 
in (|2.14p with the left state (u_,«_,6'_) replaced by {v*,u*,6*). 

Thus, from the properties of the viscous contact wave in (j2.6p and the approximate rarefac- 
tion wave in Lemma 2.3, we have the following relation between the approximate wave pattern 
(y, [/, ©)(r, y) and the exact inviscid wave pattern (V, U , G))(r, y) of the Euler equations 



|(y,C/,e)(T,y)-(y,i7,G)(T,y)|< 



c 



1 + T 



ln(l + r) + C(5'-^e"T+^. 



(2.17) 



Hence, to prove the zero dissipation limit (j2.3p on the set Sj^ defined in (j2.4p . it is sufficient to 
show the following time-asymptotic behavior of the solution to ()2.2p around the approximate 
wave profile (j2.16p . i.e.. 



lirn ||(u,n,6')(r, 

■— 5>+00 



(y,[/,G)(T,.) 



L°° 



0. 



(2.18) 



First, by (j2.7|) and (I2.15|) . the superposition wave profile (y, [/, G)(r, y) defined in (I2.16|) 
satisfies the following system 



f Vr-Uy = {), 



R 



e, 



[/2 



-e. + pc/, = K-^), + -^ + Q2, 



where P = p(y, 0) and 



P 



CD 



P«3)^ 



c/„ 






Qa ={PUy-P^' U^^ - P^^ U^^ - P^^ [7^3 ) - 



V 



yCD 



Q 



CD 



V 



yCD 
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(2.19) 



A direct calculation shows that 

+\{yy^',Qy')\\{y^'-v*,Q^'-G*,v^^-v*,Q^^-e*)\ 
M{u;;'',v^''muy^\Vy''\Uy^^v^n\ + my^\y^'muy''^yy''n\} 

'■= Qu + Qi2- 

Similarly, we have 

Q2 = o(i){|c/f ||(y^^-t;„e^^-0„y«3-7;*,e^3-r)| 
+\u^^\\{v^^ -v,,e^' -e,,v^^ -v*,Q^^ -e*)\ 
MiUy'', yy^"", e^^)||(y^^ - v., 6^1 - e,, y^=^ - v*, e^^* - e*)\ 
+|(^,^^,Kf^,G^^)||(c/,«SKfsef,t/«Mf«,ef)| (2.20) 

+l(c/f,<^ef)ll(c/f,vf^ef)|} 
+o(i){|e«M + |e«^| + |(c/,^Sif\ef,f/f,y«3,ef)p} + |Q^^I 

Here Qu and (521 represent the wave interaction terms coming from the wave patterns in the 
different family, Q12 and Q22 stand for the error terms due to the inviscid approximate rarefaction 
wave profiles, and Q is the error term defined in (12. Sp due to the viscous contact wave. 

In fact, one can estimate the interaction terms Qu and Q21 by dividing the whole domain 
^ = {(''"! y) I (''") y) £ R- X R} into three regions: 

f]_={(T,y)|2y<Ai,(l + r)}, 

^CD = {(r, y) I Ai.(l + T)<2y< A^(l + r)}, 

n+ = {{T,y)\2y>Xl{l + T)}, 

where Ai* = Ai(v^,,0*) and A3 = X3{v*,9*). Then, in each section the following estimates follow 
from (12. 6p and Lemma 12.21 

• In il_, 

• In O.CD, 

|(y«i _ ^„vfi)| = 0(i)<5Rie-2{M+|Ai.|(i+r)}^ 

|(y«3 _^*^V;iJ3)| = 0(l)5iJ3e-2{|x| + |A5|(l+r)}. 
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Inn 



+) 

Keep in mind that each individual wave strength is controlled by the total wave strength by 

([LB]) and ([L7]), that is, 

§Ri + ^^3 + §CD < (^^_ 

Hence, in summary, it follows from (j2.19p . (|2.20p and the above arguments that 

|(Qii,Q2i)| = 0(l)5e-^^l^l+(i+^», 

for some positive constant C independent of r and y. 

3 Proof of the main result 

In this section, we shall prove the main result Theorem 11.11 By virtue of the arguments in 
Section 2.3, it is sufficient to show ()2.18p besides the regularity of the solution. To this end, we 
first reformulate the problem. 

3.1 Reformulation of the problem 

Set the perturbation around the wave profile (V, U, 0)(r, y) by 

((^, i;,C){T,y) = {v,u, e){T, y) - {V, U, 9) (r, y) . 
Then, after a straightforward calculation, the perturbation {(p,ip,Q{T,y) satisfies the system 

(pr-'4'y = 0, 

D R 9) V? U'^ (3.1) 

^C. + (Pn, - PU,) = K^ - ^)y + (f - f ) - Q., 

(0,V',C)(7" = 0,y) = (0o,V'o,Co)(y), 

where the initial data ((/>o,'0Oi Co)(y) and its derivatives are sufficiently smooth away from but 
up to y = 0, and 

('Ao,V'o,Co)(y) G L\n),ct>^y G l2(r±). 

For simplicity, denote 

In order to prove (j2.18p . we easily see that it suffices to show 

12 



Proposition 3.1. There exists a positive constant 5o, such that if the wave strength 6 and the 
initial data satisfy 

S+Afo <So, 

then the problem (j3.ip admits a unique global solution {(j)^%l)^C,){t^y) satisfying 
(i) There exists a positive constant C independent of t, such that 



sup 

r>0 



\ r^"^ 1 



(a) For any tq > 0, there exists a positive constant C = C(ro), such that 



Xy,i^r,Cr)ir,-)f + 



sup II \^y, <,y, 

r>ro J TO 



+ 00 



Cyy,i^yT,CyT){T,-)f dr < C{to){Mo + 6^). 



yy, Syyi Vyri Syr 



(Hi) The jump condition of (J){t, y) at y = admits the bound 

|[0](r)|<Ce— 
where the positive constants C and c are independent of t £ (0,+oo). 
Assume that Proposition 13.11 holds, then for any tq > 0, one has 



(3.2) 



+00 



TO 



d 



y^i^yXy)^^ +1 ^11 (</>?/, V^j/, C?/)-[f^ Ijdr < +oo, 



whence, 



lim 



i^y,Cy)f=0, 



y, ^y, Sy 



which, together with Proposition 13. II and Sobolev's inequality, implies that 

lim sup||(</>,V,C)llioo <C7 lim ||((^, V', OIIIK'^y, V'y,Cj;)^ < C lim ||(<^j„ ^„ C,)^ = 0. 

T^OO ^g T— >00 a a a T^tOO 



The above inequality combined with ()3.2p gives ()2.18p . Thus, the main result Theorem 11.11 
follows from (f2T8|) and (|2T7J) . 
Denote 



iV(n,r* 
and define the solution space by 



sup {\m^,C)ir,-)f + \\{^y,^yXy){r,-)f +\\i^r,Cr){rr)f}, 

t6[t.,t*] •- J 



X[n,T*] = { 



(</',^,C) 



l/-D±\ 



,V',C)(r,y)GC([T.,T*];iJ^(R*)) 



(^T,Cr)G^°°(n,r*;L2(R±))nL2(T,,r*;Fi(R±)). 



Since the local existence of solutions to (|3.ip is proved in [7], we just state it and omit its 
proof for brevity. 
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Proposition 3.2. (Local existence) Suppose that Mq and the wave strength 5 are suitably small 
such that inf fo 0''nd voIOq are positive. Then there exists a positive time tq = tq{N{<S)^5) > 0, 
such that the Cauchy problem (j3.ip admits a unique solution {(j),'ip,Q{T,y) € X[0,ro] satisfying 



AiTo) + B{to) + F(ro) < C(AAo + 6), 



where 



A{to)= sup |||(</>,^,C)(r,.)f + ||<^,^'}+ TlK^^Ofdr, 

0<r<To "- ^ JO 



B{to)= sup \g{T)2 

0<T<TO ^ 

;.2||2 



+ 9ir)\\4>y^'}+ rg{T)h^''{Urf + \\i^)yf)dT 
' Jo V 



+ 



9{r)m 







yi 



+ 



+ \\(^)yf)dr, 



F{to)= sup \g{T)2^ 

0<T<TO ^ 

f-ra , 

+ / giry^^'^Uyrf +9{Tf\\Cyrf)dr, 

Jo 



+ \\i^)yf) + 9{r)H\\Cr\\' + \\{\f)} 

V V ■) 



p- 


_Uy^ 


= 


A- 




V . 




- V . 



with g{T) = r A 1 = iiiin{r, 1} and "d G (0, 1). Moreover, v,u,6 have the same regularity as in 
Theorem \1.1[ Thus, v,Ux,Ox have one-side limit at y = and satisfy the jump conditions 

0. 

Finally, one has the following estimate on the jump at y = 0, 

\[v]{T)\<C6e-^\ T>0 

for some positive constants C and c independent of t. 

Hence, in view of tlie local existence and the standard continuation process, we see that to 
prove Proposition 13. H it suffices to show the following (uniform) a priori estimate. 

Proposition 3.3. (A priori estimate) Suppose that the Cauchy problem ()3.ip has a solution 
{4),il),C,){T,y) E X[ti,T2]. There exists a positive constant rji, such that if 



iV(Tl,r2) + 5<7?l, 



(3.3) 



then, 



iV(Ti,T2)+ r{||./.,(r,.)^2+||(^„C,)(T,-)^f+||(V>Cyr)#(r,-)}rfr<C(iV(Ti) + 5i), (3.4) 

Jti 

where the positive constant C is independent of t. 
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3.2 Energy estimates 

In this section we will derive the a priori estimate given in Proposition 13.31 Note that under 
the a priori assumption ()3.3p . if ?y ^ 1, then if holds that 

inf {{v + 4>,e + 0{r,y)}>Co 

[ti,T2]xR 

for some positive constant Cq. First, one has the following Lemma: 

Lemma 3.4. Under the assumptions of Proposition \3.3[ there exists a constant C > 0, such 
that for any r G [ti, T2], 



Jti 

<C\\{cP,i;,C,<PynHri) + C ril + Trl\\{cP,i;,C)i;r)fdT + C6'^ 

Jti 

+ C5 r -f {l + T)-^e~^^\{cP,C)\^dydT. 

Jti Jr 

Proof: Let 

$(2) = z - 1 -Inz. 

Arguing similarly to that in |12j or J14j . one can get the following equality 



^y . ,._l^ I pCf/^i +[7R3^U.f_■ 



= Qi- Qi'>P-Q2-x, 



where 



and 



(3.5) 



/r(r,.) = i?e<.(^) + f + :^^(|), 

^i(r,y) = (p-m-(^-^)^-K^-^)^, (3.6) 



Q3 = -PU-- (<l>(^) + 7'^(^)) + (.(^), + I + Q.) {(7 - im'y) 

+M-) - ^1 - (- - -)f/ ^ + (i - -)U'^- + 2^^^ + .^^ (3.7) 
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Integration of the equality (j3.5p with respect to y and r over R='= x [ri , r] yields that 

',2 CV/-2- 



Ii{Ti,y)dy+ i {Q3 - Qii) - Q2-^)dydT. 



It is easy to observe that the jump of Hi in ()3.6p across y = vanishes, i.e., 

= [p-^]«(..0)-[P-f]^.(.,0)-.([^]-[|])|^ 

Recalling that 

$(1) = $'(1) = 0, $"(z) = z"2 > 0, 

there exists a positive constant C, such that if z is near 1, then 

C-\z-lf <<^{z)<C{z-lf. 

Thus under the a priori assumptions (j3.3p . one gets 



coy 



(3.8) 



C-^P < ^(f ) < C|.^|2, C-1|CP < <1'(q) < C|CP (3.9) 

and 

C7-i|(0,C)P < ^(^) + 7'&(:^) < C\{cP,0\'. (3.10) 

Now it follows from (13. 7p . (13. 9p . (|3.10p and Cauchy-Schwarz's inequality that 

+(|(y/^3^c/f,ef)|Me,«,3|) + |g2l|(^2^^2), 



(3.11) 



i = 1,3, 



By the properties of the viscous contact wave, one can obtain 

r / (le^^p, le^jf |)(02 + c')dydT < C5 r Ul+ r)"ie-iTF|(<^, 0\'dydr, 
while by the properties of the approximate rarefaction wave in Lemma 12. 2| we have that for 

1 JR. 

Jti 

<C f\l + Tr'\\{^,0\\my,Cy)\\dT 

</. / \\{cPy,CyWdT + C^ ril + ryWicPXWdr, 
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where and in the sequel /i is a small positive constant to be determined and C^ is some positive 
constant depending on fi. 

Now, it remains to estimate the terms Qiip, Q20 on the right-hand side of (j3.8p and the 
term \Q2\{4>'^ + C^) on the right-hand side of (|3.11|) . For simplicity, we only estimate Q2^- By 
(122nD, we find that 

[ j^\Q2'^^\dydT <C j^^ \\Q\\L^\\Q2\\LldT 

< cj^ IICII^IIC,II^(||Q2i||li + \\Q22\\li + \\Q''''\\Li)dT 

Jti 

<fi r wcyfdT + c^ 6-e r iiciii(i + r)-idT 

< ^1 r WCyfdT + C^ r llCf (1 + r)-ldT + C^ <53. 



'n 



Similarly, one can control the term Qiip and |<52|(</'^ + C^)- 

Thus, substituting all the above estimates into (j3.8p and choosing /i in the front of the 

integral / IKV'yi Cj/)ll dr small enough, so that the integral can be absorbed by the left-hand 

side of (|3.8p . one concludes 



• i^X){rr)f+ n\\{^y,Cy){r,-)f + \\^J{Uy'^'\Uy'''n{^,C){r,-)f}dr 
<C||(0,V,C)(Ti,-)f + C / il + T)-l\\i<p,i;X)fdT + C6-^ (3.12) 

+C/X r\\cl,y{T,.)fdT + C5 r Ul + T)-^e-'^UX)?dydT. 

Jti Jti Jr 

Next, we estimate ||(/>y|p. Denote v = y. From the system ()3.iP o. one has 

(^)r -A-{P- P)y - Ql = 0. 



V 

Vy 



Multiplying the above equation by -y- and noticing that 



V V V V V V 



one obtains 



^ ,Vy 2 i'^y\ f 1 ^t\ ^^/^y\2 



j^r-i^^ +u^ +— (^) 

2 V V J ^ \ V J y V V 

V V \ V V V V I V 
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y 



Integrating the above equality with respect to y and r over H^ x [ti,t] and using Cauchy- 
Schwarz's inequality, we infer that 



< 






V V ' 



V'— 

V 



{T)dT + 



ff 

Jti Jr 



Re ,u 



y\2 



2v V 



dydr 



V V 



+C 



Tl JR 



RCy 



{Ti,y)dy+ / f \ipy{ 
Jti Jr. 



"'7/ ^1i 



V 



Qi 



)\dydT 



dydr, 



where the jump across y = can be bounded as follows. 



V 



{T)dT 



^(r,0) 



"y 

V 



'y 
V 



{T)dT 



^P{T,0)[p]{T)dT 



R / V(t,O)0(t,O) 



{T)dT 



-rT :»':'°"''"-°'.M(.)d. 



v{t,0+)v{t,0-) 



<c 



<6 



Loo(r)|H|(Ti)e-^(^^^i)(ir < C5 
\'^dT + 6 sup llV'f (r) + C(5. 



|ig-C(r-ri)^^ 



Using the equality 



t6 Ti,r2 



Vy _ Vy Vy 



by Vy^ 



V V V vV 



we see that 



C-\\<Py\' - \Vy<P\') < C^f < Ci\cPy\' + \VyCl>\'). 



From the definition of Qi in (j2.19p it follows that 

y" WQifdT < C j' (lIQnf + \\Qi2\?)dT 

Therefore, substituting all the above estimates into (13.13^ . we conclude that 

ll'/'.(r,-)^'+ ril</'.^'dr<C||(0,V^,0,)#(Ti) + C||(,^,^)(T,-)f 

Jti ^ 

+C5 r I {l + T)-'e-'^\{<P,C)\'dydT + C [ \\{^yXy)fdT 

Jti Jr Jti 

+C r{l + Trl\U,ij,Q)fdT + C5'^. 



Tl 



(3.13) 



(3.14) 



Multiplying the inequality (j3.12p by a large constant Ci > 0, and summing the resulting 
inequality with (|3.14p . we obtain Lemma |3.4[ This completes the proof. D 

Next, we derive the higher order estimates, which are summarized in the following Lemma: 



Lemma 3.5. Under the assumptions of Proposition [3731 it holds that 

'■T2 



N{n,r2) + / ' {\\\/(ufjjf)i^X)f + Uyf +\\i^y,Cynl +\\i^yr,Cyr)f }dr 

<CN{n) + C r\l + T)-l\\{<t^,i;X)fdT + C5\+C5 r l{l + T)-^e-'^\{4>X)?dydT. 

Jti Jti Jk 

Proof: Multiplying the equation (jS.iP o by —'4'yy, one gets 

(f) -i^ri^y)y + ^ = {iP-P)y+"f.^y-{Uyil-l,)), + Qi}^yy 
Integration of the above equation with respect to y and r over R x [ri, r] gives 
/ ^iT,y)dy+ r S ty]LdydT = -f ^(ri,y)dy- f [^r^y\{T)dT 

Jn ^ Jti JR ^ JR ^ Jti 

^ f^{iP-P)y + 5^. - {Uy{- - y))y + Ql}^yydydT =: ^ J,. 

We have to estimate Jj. First, the jump J2 can be bounded as follows. 



J2 



[ipri'y] (r)dr = - 
V'r(r,0) [Uy] {T)dT 



V'r(T,0)[V.i;](r)dr 



,u,, 



^,(r,0) {^-p)v {T)dT 
L V J 



,u. 



MrM— - P){r,Q)[y\{r)dT 



< C I \Wt\\l 

< C5 



1 1 

yr 



2 lh/;,„-[f"2 






s/ii iiv^rai 



In view of (|3.iP a and ()3.3p . one has 



< C7 

< C7 



^yy^ +II('/'J/, ^V. Cs/)^ +ll(f^yj;, ^y, f^y, ej,) 
'jyyl+Wy,i^yXy)l+^)- 



+ 



IIQill) 



Substituting (I3.17P into (I3.16|) . we obtain 

IJ2I < C8 r (\\^l;yy^+\\{cPy,^PyXy)^+sY'UyT^^ ^ 1 1 ^^y 1 1 ^ 1 1 V-,, ^ + 1^ e- ^ ^^ ' ^^^T 



< 



^ 



TO,V'yr)-[f^ dT + C^5 



y>V'y,Cj/)-[f^ dT + C^6. 



(3.15) 



(3.16) 



(3.17) 



(3.18) 
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On the other hand, J3 can be estimates as fohows. 

<C J^^ £[\{cPy,Cy)\ + \i<l>,C)my,Vy,ey,Uyy)\ 

+ \i(py,Vy)\\{i^y,Uy,Uy^)\ + \Ql\'j\ijyy\dydT (3.19) 

<l^ r UyyfdT + C^ r\\{^y,^yXy)fdT + C^ /V + t)"^ IK*/-, V', C)f ^r 

+ C^6 + C^5 r -f {1 + r)-ie-?^ !(</., Ol^dydr. 
jti Jn 

Substituting (j3.18p and (|3.19|) into (j3.15p and choosing /i suitably small in the front of the 
integral J^ llV'rall^ dr, we deduce that 

Uy\\\r) + r Uyyi^ dr < CUy\\\n) + Cf, r Uyrf dT 
J Tl J Tl 

+C^ r{l + rrl\\{^,C)fdT + C^6 + C^ r \\i^y,^yXy)f dr (3.20) 

J Tl J Tl 

+C^5 r -f {l + Ty^e-'^\{<P,C)\^dydT. 

Jti Jb. 

Multiplication of the equation ()3.iP o with —Cyy yields that 



3 "^"^^ ^yy 
y \ I r /- \ \ ,,^yy 



All) -^^^^^y\^-- 

■J — 1 \ 2 I 7— 1 ^ V 



(p-y - PUy) + .^ - .(e,(i - i))^ - (f - f ) + Q2}c, 



yy 

± 



the same arguments as those used for llV'yll^ (''") in (|3.2Up . we obtain 



Integrating the above equality with respect to y and r over R x [Ti,r], and employing almost 



(3.21) 



\\Cy\\\r)+ \\CyyiUT<C\\Cy\\\n)+C6-^+C (l + T)"a||(0,C)fdT 

J Tl J Tl 

+C ^\\{cPy,^|;y,CyWdT + Ci5f T / (1 + ^)"'e-^|((/., Ol'dydr, 

Jti Jti Jr 

where we have used the following jump estimate across y = 

-^ r [^-^J/] (^)^^ = "A rCr(r,0) [Cy] {T)dr 

I ^ J Tl I ^ J Tl 

= -- ^ rC.(T,O)[0,](r)dr = — ^ rCr(r,0)^(T,0)M(T)dT 



< c / ||Ct||l-(i + IICj;llL-)H(n)e-^("^"^)dT 

<C6 [ ||C.# ||C..# (1 + IIC.II^IIC..# )e-^(--^^)dr 



Tl 



20 



and the estimate 



\\Cr\\ < C(\\Cyy^+\\{^y,^y,(yn+my,eyy,eyVy,U^){cl,,C)\\ +\\Q2\ 



<C[\\Cyy^+\\i^y,i^y,Cyn+S 

It follows from ([HTTl) and (|3:22D that 

Cr){r,-)fdr 



(3.22) 



r2 



Tl 



<c 



yy^QyJu "t + 



n 



V's/,Cy)# d-r 



y, Yy^sy 



+ / ||([/„G,„G,y„C/2)(^,^)||2rf^+ / ||Q2||2rf^ 



<c 



yyi "^yyj 



Xyy)rdT+c \\(cPy,i;y,Cy)rdT+ / ii+T)-'\\icpx)rdT+c6-^. 



(3.23) 
Now we turn to control sup IKV't) Cr)"[f • First, applying the operator dr to the equation 

TG[Ti,r2] 

([3TT|) 9. we get 

Multiplication of the above equation by tpr gives 



2 /^ V 



2 



K(^-p).-^.(^-^).} 



y 



-Ipyr-^ + V'j/r^^^r + ^yT(-r7)r + tpyriP " -P)r " IprQlr- 



^yr^2 



If we integrate the above equality with respect to y and r over R^ x [n, r], we find that 



i4(r,y)dy+ ri^dydr 



i'l 



in,y)dy 



R 



,Uy 

V 



.[/„ 



M'^-P)r-Mf-P) 



{T)dT 



(3.24) 



+ 



71 

1 JR •- 



i^yr-^ + i^yr-^Vr + '^yr{-^)T + V'j/r(P " -P)r " i^rQlrjdydl 



V 



-S/-^2 



where the jump across y = in fact vanishes, i.e., 



U^ 



M^-P)r-M^-P)ri- 



V 



,Uy 

V 



Uy 

V 



mir){^-p)^iT,0-)+Mr,0+) i^-p)^ ir)-[A]ir){-^-P)^ir,0) 



■ U„ 



V 



(3.25) 



]r{r){^-p)^{T,0-)+Mr,0+) 



n„ 



U„ 



.(T)-[V]r(r)(-f-P)>,0) 



V 
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Now we apply dr to the equation dSJjls to deduce that 



R 



■Ctt = J^( 



,e,, 



-^ — 1 ^ V ^y^ ^ V ^y^ 

Multiplying the above equation by ^t-, one has 

Sl/T 



+Ke-p)}, 



,a 



y 



^yl-f -^)}^-Q2r. 



R ,C 



e„ 



7- 



1(1)^ + "^ = {"<-(?). -"<'(?).} 



kiyr 



e„ 



,u„ 



+ vQyr-^ + Kyr^Vr + KyTi-rj-)T + Ct%t(— " p) 

V v^ V V 

+CrUy{^ - p)r - CrUyri^ - P) - CrUy{^ - P\ - CrQ2r. 

Integrating the above equality with respect to y and r over R='= x [ti , t] , we deduce that 






,e. 



V 



^Cyr—^ + Kyr^Vr + i^CyT(-r7-)7 
f u^ V 



(3.26) 



V V 



P)r - CrUyri^ " P) " CrUy{^ - P)r " CrQ2r]dydT, 



where the jump in fact vanishes 



WC.](r)(^),(r,0-)+Kr(r,0+) 



V 



'^[C]r(r)(^)^(r,0-)+K.(T,0+) 
0. 






e, 



(r)-z^[Cr](T)(^)^(r,0) 



F 



e. 



(3.27) 



W-WC].(T)(^),(r,0) 



y 



Hence, taking into account p.25p and (I3.27p . we get from (I3.24p and (I3.26P that 



.Cr)f{t) + 



,r,Cj/r)r dr < C||(V.,Cr)r (n) + C / ||(Vr,Cr)r ^T 



+c ni+T)-i\\{cp,c)fdT+c6+c r 

+C8 r i {l + T)'^e-^-^\{(t)X)?dydT. 
Jti Jr 



,i^v,Cv)\\ dr 



y, YyiSy 



(3.28) 



Combing the estimates (I3.20|) . (I3.2ip . (I3.23p . (I3.28P and Lemma [3741 together, we obtain 
Lemma |3.5|, and the proof is completed. D 

It remains to control the term 

5 r Ul + Tr^e~'^\{(t^,C)?dydT, 
Jti Jr 
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which conies from the viscous contact wave. We shall use the estimate on the heat kernel in [12 
to get the desired estimates. 

Lemma 3.6. Suppose that Z{t, y) satisfies 

ZGL°°(0,r;L2(R±)), Zy e L'^{0,T;L^{K^)), Z^ e L^{0,T;H-^{'R^)), 

then 



■iPy 



l + r)"^Z^e ^+r dydT 
T-i Jn. 



<Cp!^\\Z{n,y)f + £\\hyfdT + l\Zr,Zgj)H,^H-Hn±)dT^ 



(3.29) 



'n 
where 



1 fy Hv'^ 
(l + r)"2 / e ^+^dr] (3.30) 

Jo 
and (3 > is the constant to be determined. 



9l3{T-,y) 

'0 



Remark 3.7. Lemma \3.6\ can be shown using arguments similar to those in US^ . and hence its 
proof will be omitted here for simplicity. Note that the domain considered here consists of two 
half lines R^ , and hence the jump across y = should be treated. In view of this, the functional 
gi3 should be chosen in (j3.30p . so that gp is continuous at y = 0. Furthermore, it holds that 
91s{t,0) = 0. 

Lemma 3.8. Under the assumptions of Proposition [XM it holds that 

_£oy_ 
V ^^l('/',^,C)Prf?/dr<C<5 + C7||(0,V,C)(Ti,-)f + C||(</.,^,C)(r,-)f 

+ Cr\\icl)y,i;yXy)fdT + Cril + T)-l\\{ct>,ij)\\''dT. 



Tl 



Proof: Prom the equation p.ip n and the fact p — P = — — - one gets 

v., + (^^— ^), = (^ - ^)y - Qi. (3.31) 



Let 



ry ^ 
Ga{r,y) = {l + Ty^ / e^^dij, 
Jo 



where a is a positive constant to be determined. Multiplying the equation (j3.3ip by Ga{RC- 
we find that 

S j,-<°"'»^^ + ^' <''»+« (3.3,, 

-Ga{RC - P(I))A + Ga{RC - P(l)){^ - -^)y " Ga{RC " P^)Ql. 
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Noticing that 

- G^{RC - P(^)i,r = -{GaiRC - P(^)i')^ + {Ga)r{RC " ^0)^ + G^Vl^C " ^0)r (3.33) 

and 

{RC - P(p)r = RQr - PT(t> - P(t>T 

= -7PVy + (7 - 1){ - (P - ^)(f/. + ^.) + (f - -^) + Kf - -^)y - Q2} - P.0, 
if we insert p.34p into (j3.33p and use the equality 

we get from ()3.32p that 

ay 

2(l + r) {^^'^ " ^"^^^ + ^^"^'1 " ^^"''(^^ - ^"^^^Ir + ^2, + Q4, (3.35) 

where 

2v 2 V V V V 

and 



= ^°^^2.2 ^'^^ V . + 0.) - iGa)r{RC - P^)^P + (G„(i?C - P0)),(^ - ^ 

+(7 - 1)G„V' (p - P)iUy + ^,) - (-1 - -^) + Q2 



yK^ - ^) + G,(i?C - i^'/')Ql + Ga^Pr^ - jPy^ 



+(7 - l)z/(G„V)y(^ " "/) + ^"^^^ " ^'^^^i + G^V'^r'/' - 7^J 
Integrating (j3.35p over R^ x [ri , r] , one infers that 



/ / ^^{(^C-^'/')' + ^'}rfy^r = / {G,t-(iiC-^0)^}(r,y)dy 
-/ {GaKi?C-^0)^}(n,y)rfy+ [ [H2]{T)dT + f i Q^dydr. 



(3.36) 



Here we only analyze the jump term [H2] across y = 0, the other terms in (|3.36p can be estimated 
similarly to those in [12] or [14] . Recalling that Ga{T, y) is continuous at y = and Ga{T, 0) = 0, 
we easily see that 

[H2]{r) = 0. 
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Thus, from p.36p one gets 



n JK 



_ay_ 

e 1+^ 

1 + T 



{{RC - PcPf + i>^]dydT <C5 + C||(</), V,C)(ti, 



+C||(0,V,C)(T,-)f +C/ (l + T)-6||(0,V^,C)(T,-)fdr (3.37) 

ay 

+C r\\{<Py,^y,Cy){r,-)fdT + C6 r -f ^^|(0,C)|2dydr. 

In order to get the desired esthnate in Lemma 13.81 we will use Lemma 13.61 to derive another 
similar estimate from the energy equation (j3.iP n. To this end, we set 



Z 



R 



7-1 



C + Pcj) 



in Lemma [3.61 Thus we only need to compute the last term in ()3.29p . Prom the energy equation 
()3.ip n. we have 



Zr = Pr(t> - {P- P)Uy + 1^{ 



V V 



+ 



' V V ' 



Q2 



whence 



Jti 

r / {Pr4> -{p- P)Uy)ZgjdydT - r -f {p- P)^yZg}dyd 

Jti Jn Jti Jr 

^ .(^ - %Zg}\ {r)dr - fl ui^ " %{^Zgl)ydyd, 
f i^- -f)ZgldydT - j j Q2Zg}dydT =: ^K,. 



Here the jump term K^ can be estimated as follows, recalling 5/3 (r, 0) = 0. 



K3 



n 



K 



^y ®y 
V V 



)Zg$ {T)dT = u g'p{T,0){ 



n 



'^-^-^){r,0)[Z]{r)dr^O, 



while the terms Ki {i = 1,4,5,6) can be directly dealt with in the same manner as in [12] or 
. To bound the term K2, we make use of the mass equation (j3.ip ^ to write K2 in the form 



{P - P)^yZg 



a 



lP^-{l-l)Z^ _-lPZg^ 



Zgp(pT 



(7 - l)Z'gl 



V '' 2v 

-tPZ<P^g} - 2(7 - 1)0^2^2 ^p^^2 _ 2(^ _ 1)^2^ 



2v 



V 



-9l3i9l3)T 



^jPZ£-2(j-l)Z^ 2 /2(7-l)5^ tF^n 19jrz 



2^2 



2v 



2v 
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where all terms on the right-side hand of the above identity can be directly bounded in the same 
way as in [12] or |14j . Therefore, we have bounded i^2- 
Taking /? = y-, one can get from Lemma 13.61 that 

f I 2^z^dydT<C5 + C\\icP,^PX){ri,-)f + C\\{^,i^,C)ir,-)f 
+C r\\i<Py,,Py,Cy)fdT + C ril + T)-l\\i<p,i;)fdT (3.38) 

+C{6 + m) r / (1 + T)-ie-^|(0, Ol^dydT. 

Jti Jr 

Now, taking a = co in (I3.37P and choosing 5 and r/i suitably small, we combine (I3.37P with 
p.38p to obtain the desired estimate in Lemma 13.81 D 

By Lemmas 13.51 and 13. 8| we conclude 

N{n,T2) + r {Uyf +\myXy)€ +\myr,Cyr)f }dT 



<CN{n) + C [ {l + Trl\\{^,^,0fdT + C6l 
Jo 

An application of Gronwall's inequality to the above inequality gives the estimate ()3.4p in Propo- 
sition 13.31 This completes the proof of Proposition 13. 1[ 
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